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Mixtures of relativistic gases in gravitational fields: 
combined Chapman-Enskog and Grad method and the Onsager relations 

Valdemar Morattc0 and Gilberto M. Kremei0 
Departamento de Ftsica, Universidade Federal do Parana, 81531-980 Curitiba, Brazil 

In this work we study a r-species mixture of gases within the relativistic kinetic theory point 
of view. We use the relativistic covariant Boltzmann equation and incorporate the Schwarzschild 
metric. The method of solution of the Boltzmann equation is a combination of the Chapman- 
Enskog and Grad representations. The thermodynamic four-fluxes are expressed as functions of the 
thermodynamic forces so that the generalized expressions for the Navier-Stokes, Pick and Fourier 
laws are obtained. The constitutive equations for the diffusion and heat four-fluxes of the mixture 
are functions of thermal and diffusion generalized forces which depend on the acceleration and the 
gravitational potential gradient. While this dependence is of relativistic nature for the thermal force, 
this is not the case for the diffusion forces. We show also that the matrix of diffusion coefficients 
is symmetric, implying that the thermal-diffusion equals the diffusion-thermal effect, proving the 
Onsager reciprocity relations. The entropy four-flow of the mixture is also expressed in terms of the 
thermal and diffusion generalized forces, so that its dependence on the acceleration and gravitational 
potential gradient is also determined. 

PACS numbers: 51.10.-|-y, 05.20.Dd, 47.75.+f 


I. INTRODUCTION 


The relativistic kinetic theory of gases is a subject that 
began in 1911 when Jiittner [l| proposed a relativistic 
version of the velocity distribution function which cor¬ 
responds to the Maxwellian distribution function in the 
non-relativistic limiting case. Later, several studies have 
been made, but for brevity’s sake we mention the books 
ii where several applications of the relativistic kinetic 
theory of gases are discussed. 

This work represents a continuation of the study of 
the properties of relativistic gases using the Boltzmann 
equation in gravitational fields, this sub-area has noWet 
been studied in depth. Here we quote some works 
on this topic which have been recently published. 

The method used in this paper to solve the covariant 
Boltzmann equation is a combination of the Chapman- 
Enskog and Grad methods [^, [l^ . It consists essentially 
in doing an expansion of the distribution function for 
each species which is the solution of the Boltzmann equa¬ 
tion up to first correction. Such a procedure is like in the 
Chapman-Enskog method. Then we impose that such 
an expansion must be compatible with the solution of 
the Boltzmann equation given by the method of Grad 
11| . In order to keep the linear regime we truncate the 
Grad distribution function for each species up to linear 
terms of the non-equilibrium pressure, pressure deviator 
tensor, diffusion and heat four-fluxes. Therefore, we ob¬ 
tain a linearized Boltzmann equation that is written in 
terms of the local thermodynamic variables and fluxes: 
diffusion, heat, non-equilibrium pressure and pressure de¬ 
viator tensor. The next step is to get from that linearized 
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Boltzmann equation a set of linear algebraic system for 
the fluxes. We generate one equation for each thermo¬ 
dynamic flux through the multiplication of a dynami¬ 
cal function of the particles by the linearized Boltzmann 
equation and then the integration over the momentum 
space. Hence, we find the constitutive equations for the 
fluxes in terms of gradients of the local thermodynamic 
variables and of a gravitational potential that arises from 
the Schwarzschild metric. The laws of Navier-Stokes for 
the non-equilibrium pressure and pressure deviator ten¬ 
sor are obtained as well as the generalized Fourier and 
Fick laws for the heat and diffusion four-fluxes. 

It will be shown that there appears a generalized diffu¬ 
sion force that has, not only dependence on the concen¬ 
tration and pressure gradients but also on a contribution 
of the four-acceleration and the gravitational potential 
gradient. The contributions of four-acceleration and po¬ 
tential gradient also appear as combined forces for the 
heat flux and the y w ere analyzed separately by Eckart 
[l^ and Tolman |l3l . M- In the case of Eckart, for a 
relativistic gas in equilibrium and in the absence of grav¬ 
itational helds, the temperature gradient is counterbal¬ 
anced by an acceleration. On the other hand, in the case 
of Tolman for a relativistic gas in equilibrium and in the 
absence of an acceleration, the temperature gradient is 
counterbalanced by a gravitational potential gradient. 

In order to show that the Onsager reciprocity rela¬ 
tions hold we manipulate the constitutive equations for 
the heat and diffusion fluxes. The demonstration is gen¬ 
eral in the sense that the interaction of the particles are 
supposed to maintain the microscopic reversibility prin¬ 
ciple. 

The structure of this paper is as follows. We define 
the problem in section H and establish the Boltzmann 
equation and the definitions for both, the thermodynamic 
variables and fluxes. In section HI, we use a method of 
solution of the Boltzmann equation that is a combina- 
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tion of the Grad and Chapman-Enskog ones, the solu¬ 
tion is truncated up to first order so we obtain linear 
expressions. Such a process will lead us to an algebraic 
system of equations for the thermodynamic fluxes that 
when it is properly solved, expresses the thermodynamic 
coefficients for an arbitrary inter-molecular interaction. 
In section IV, we show that the Onsager reciprocity re¬ 
lations hold for an arbitrary inter-molecular interaction. 
Furthermore, we show that the laws of Fourier and Pick 
are expressed in terms of generalized thermal and diffu¬ 
sion forces in the presence of gravitational fields. To give 
a more general representation, we show in section V that 
the entropy four-flow of the mixture is a function of the 
generalized thermal and diffusion forces. Section VI is 
devoted to the calculation of the constitutive equation 
for a relativistic Newtonian fluid, i.e., the Navier-Stokes 
law. Ultimately, in section VII, we discuss the obtained 
results. 


II. BACKGROUND 


In this section we will define the problem of a r-species 
non-reacting mixture in a Riemannian space with metric 
tensor . The particles are supposed not have inter¬ 
nal degrees of freedom. Each of these particles of the 
constituent a = have mass nia and are character¬ 

ized by the space-time coordinates = (ct, x) and the 
momentum = (PaiPa)- The mass-shell condition, i.e. 
gfivPaPa = rn^c^ implies the following relationships for 
the contravariant and covariant temporal components, 

iPa = (Pa O - 90iPa)/900, __ 

\ Pa 0 = y 9oomlc^ + [goigoj - 9oo9ij) Plpi, 


respectively. The analysis is developed within the tenets 
of the general relativity, we adopt the Schwarzschild met¬ 
ric g^'^ in which the line element reads [l5l| : 


2 / 2GM \ . g, 2_1_ 2 

-R^ [de^ + sin^ Ody^) , (2) 

in terms of the spherical coordinates = x°}. 

Above, M is the total mass of the spherical source and 
G is the gravitational constant. Here we shall use the 
isotropic Schwarzschild metric, which reads 


goir) 


ds^ = go{f) {dx^Y ~ gi{i')5ijdx'‘dx\ 




gi{r) = 1 + 


GM 


2c^r ) 


( 3 ) 

( 4 ) 


Along the calculation we will use a relativistic param- 
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eter , where c is the speed of light, k the Boltz¬ 

mann constant and T the local temperature, assumed as 
an invariant. This parameter is convenient because it tell 


us how relativistic is the system, for example, Ca ^ I cor¬ 
responds to a non-relativistic limit. On the other hand, 
Ca <C I belongs to an ultra-relativistic limit. 

The most fundamental equation in the kinetic theory 
is the the Boltzmann equation; such an equation can be 
obtained with two hypothesis as a basis. The first one 
is that particles collide elastically and only collisions of 
pairs are taken into account. The second one implies a 
description of the system with one-particle distribution 
function, that is equivalent to think that collisions rep¬ 
resent a process that do not depend of what occurred in 
the past with the particles. This last hypothesis is also 
known in the literature as molecular chaos hypothesis. 
In our case, the Boltzmann equation reads 


P 


a f^i'PaPa 


= X] /“ fafb)Fba(Tab 


( 5 ) 


for the a-species. Here the Latin subindex denotes the 
species, note that we have one equation with the same 
structure of ([5]) for each component of the mixture a = 
1,..., r. The distribution function fa {x^,p^) has a statis¬ 
tical meaning, indeed the quantity fa {x^,pf) df'xdf’pa at 
time t, is the number of particles of the constituent a in 
the volume element between x, -x.+d^x and p^j, Pa + d^Pa- 
In equation ([5|) also appear the Christoffel symbols F^j, 

and the invariant flux Fba = {jJaPb^iY ~ 
which plays the role of the relative velocity of the non- 
relativistic Boltzmann equation. We have also the invari¬ 
ant differential elastic cross-section CTabdU for collisions 
of species a and 6, where dVL is the corresponding solid 
angle element. Integrals are made with the invariant dif- 

_ ,3 

ferential element being = det[g'''^]. In 

equation ([5]), quantities denoted with a prime are evalu¬ 
ated with the momentum of the particles after a binary 
collision occurs, i.e., /' = /(x, p^,t) and so on. The bi¬ 
nary collision is characterized by the energy-momentum 
conservation law Pa+Pb = P'a + p'b- 

Without solving Boltzmann’s equation we can obtain 
two important results. The first one, arises from the H- 
theorem and the definition of the thermodynamic vari¬ 
ables. A situation of local equilibrium means that the 
entropy four-flow production (see Eq. (I^ l vanishes at 
equilibrium. The solution of the collisional term of the 
Boltzmann equation - when it is equal to zero - is the 
well-known local equilibrium distribution function, which 
reads 


/r = 


A'KkTmlcK2 (Co) 


exp 


kT 


( 6 ) 


Here rio is the local number of particles of species a, the 
modified Bessel function of second kind is represented by 


Kn{0 


(CY r(i/ 2 ) 

\2j T{n + l/2) 


/ OO 

e-CV (y2 _ 1/2 

( 7 ) 
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and Ufi - with - is the hydrodynamical four- 

velocity. The set of local hydrodynamic variables that 
describes the local equilibrium is {rii,n^, For 

the calculations that will be performed, it is convenient 
to evaluate Eq. (jb)) in a co-moving frame, that is = 
(c/v^:0) yielding 

r(o) ng / + gi\pa\A 

AiTkTmlcK^ iCa) kT j' 

( 8 ) 

The second important result that arises from the Boltz¬ 
mann equation is the obtention of the balance equations, 
for this purpose we proceed as follows. We multiply the 
Boltzmann equation (O by the collisional invariants, that 
is, microscopic dynamical quantities that are conserved 
between collisions, i.e., 'bg -I- 'I'b = -I- 'I'J, and integrate 

the resulting equation over \/—g . The collisional 

invariants 4'g can take the value of the mass and the 
energy-momentum of the colliding particles. To obtain 
the particle four-flow balance equation for the a-species 
we take 'kg = c and integrate, this process leads to the 
conservation law 

= 0 . ( 9 ) 

Here the semicolon denotes a covariant derivative and we 

have defined 

= ( 10 ) 

J PaO 

as the particle four-flow of species a. We now introduce 
a general decomposition of in terms of the hydrody¬ 
namic four-velocity as 

= Uall ^where ^ (11) 

denotes the partial particle number density. The quantity 
Jg is a space-like vector defined as 

j;) = AOc /pC/g^, (12) 

J PaO 

and holds the property JgH^ = 0. Above, we have intro¬ 
duced the projector 

(13) 

that has the property = 0. Equation (Ha is the 

corresponding diffusion four-flux of species a of the mix¬ 
ture and by taking the sum of (HU over all the compo¬ 
nents we easily note that 

r r r 

a—1 a—1 0=1 

( 14 ) 

where the last equation implies that there exist only 
(r — 1) partial diffusion fluxes that are linearly indepen¬ 
dent for a mixture of r constituents. 


On the other hand, to obtain the balance equation for 
the energy-momentum of the a-species defined by 

Tr = C f ( 15 ) 

J PaO 

we multiply the Boltzmann equation ((S]) by the collisional 
invariant 4'g = cp^ and integrate the resulting equation 

_ ,3 

over ^-g-^- This process yields 

Tg^C = PH, (16) 

where the production term is given by 

PH = j2^[ {pT-P^a)fahnaaatdnV^g^V^g^. 

^ J PbO PaO 

( 17 ) 

Note that this equation does not represent a conservation 
law, but if we sum Eq. dni) over all species we obtain 

r 

rr = E^«=H (18) 

a—1 

that represents a conservation equation for the energy- 
momentum tensor of the mixture By 

following the decomposition of Eckart (see e.g. [2ll - [^ l. 

the energy-momentum tensor of the a-species can be 
written as 

J.,. = + hgJD 

+ (qO + hgjK) - (Pg + t^g) , (19) 

where several definitions are to be made. Eirst we can 
list the local equilibrium quantities: energy per particle 
eg, hydrostatic pressure Pg and the enthalpy per particle 
hg = Cg -I- Pg/ng. Next, the non-equilibrium quantities 
are: dynamical pressure rajg, heat four-flux qf) and pres¬ 
sure deviator tensor . They are given in terms of the 
following projections of the energy-momentum tensor of 
the a-species: 

qO + hgJO = eg = ^U^TrU,,{2Q) 

= ( 

Pg + -(22) 

The corresponding quantities for the mixture are: 

r r r 

e = Ev®“’ P = EPa> ro = E^“’(23) 

a—1 0=1 0=1 

h = E-ha, =Epi^'^^ (24) 

r 

= E (9a + haj;(). (25) 

0 = 1 


A;(AC--A^"Agqrg"",(21) 













4 


so that the energy-momentum tensor of the mixture is 
written as 

-{p + w)A^^’'+ (26) 

Note that the heat four-flux Eq. (1251) of the mixture 
has two contributions, this is in accordance with the Lin¬ 
ear Irreversible Thermodynamics 0 , one term is related 
with the partial heat flux and another with the transport 
of energy driven by diffusion. 

Another quantity which is important in the analysis of 
mixtures of relativistic gases is the entropy four-flow of 
the mixture S^, defined by 

=-A:c^ /'p^/aln(ba/a)v^—(27) 

Pao 


theory are due to Marie and Anderson and Witting (see 
e.g. Si). For the Marie model one should take the 
Eckart decomposition, while for the Anderson and Wit¬ 
ting model the Landau-Lifshitz decomposition should be 
used. 

The main problem in the kinetic theory is to find a 
solution of the Boltzmann equation ([5]), because as we 
have seen, all the above definitions can be evaluated by 
integrating functions that involve fa The equi¬ 

librium quantities can be evaluated with the local equi¬ 
librium distribution function ([5]) and read: 


Sa = kiln 


XflaC ^ J , 

Pa — ^akT^ ha — nflaC Ga^ 

4TrmlckTK2{Ca) 


Haba 


+ CaGa - U . 


(32) 

(33) 

(34) 


where ha is a constant which has inverse units of fa- Its 
balance equation is obtained through the multiplication 
of the Boltzmann equation ([S]) by —fccln(ba/a), the sub¬ 
sequent integration over the sum over all 

species, yielding 


= a > 0, (28) 


, r r 

= ?EE 


/a/&ln 


f' f' f f' f 

J aJ b I J aJ b 


a—1 6=1 

xFbaO-abdfly/^ 


fafb \fafb 


- 1 


d^Pb 

PbO 




d^Pa 

PaO 


(29) 


The quantity cr is the entropy four-flow production of the 
mixture, which is always positive semi-definite, thanks to 
the relationship (a;—1) Inx > 0 valid Va; > 0. The entropy 
four-flow of the mixture is decomposed according to: 


= nsC/^-b$'", s 


-S^Ua, 


= A;)5'^(30) 


where the quantity s is identified as the entropy per par¬ 
ticle of the mixture and <I>^ its entropy flux. The entropy 
per particle of species a is given by 


kUa 


Sa= - 


cn„ 


plffaln{hafa)V^ 


d^Pa 

PaO 


( 31 ) 


so that we have ns = ^aSa. 

In the kinetic theory of relativistic gases there exist 
two decompositions that are often used: the Eckart and 
the Landau-Lifshitz (see e.g. ii). The difference be¬ 
tween the decompositions is that the heat flux appears in 
the particle four-flow but not in the energy-momentum 
tensor in the Landau-Lifshitz decomposition, contrary to 
the Eckart one. One can take both decompositions for 
the determination of the constitutive equations and the 
results are the same. However, there are situations where 
one should apply only one of the decompositions, which 
is in the case of using BGK models of the Boltzmann col¬ 
lision operator. The model equations of the Boltzmann 
equation normally considered in the relativistic kinetic 


The chemical potential of species a is introduced through 
the Gibbs function per particle, namely, pa = Gq — Tsa -b 
Pa/'^a and by taking into account the above expressions 
we get 


Pa = kT In 


engba 

A-KmlckT K2{C,a)' 


(35) 


In next sections, we will use a method that allow to 
obtain expressions for the diffusion fluxes Jg, heat flux 
non-equilibrium pressure tu, pressure deviator tensor 
pGG and entropy flux Furthermore, we will show the 
dependence of JJ), and q^ in terms of the gravitational 
potential and demonstrate the validity of the Onsager 
reciprocity relations. 


III. COMBINED CHAPMAN-ENSKOG AND 
GRAD METHOD 

In this section we will use a method to extract ther¬ 
modynamic information from the Boltzmann equation 
[§, 0 that combines the features of the Ghapman- 
Enskog 0 and Grad’s moments one [ll|. This method 
has mainly two advantages, the first one is that we do 
not need a solution of the integro-differential Boltzmann 
equation as in the Chapman-Enskog method. The sec¬ 
ond is that we do not need the field equations for the 
moments as in the Grad method. 

First we are going to describe how does the moment 
Grad method is constructed. The central idea is to ex¬ 
pand fa (a;^, pO) around the local equilibrium distribution 
function in a series of an ortho-normal set. In this case, 
we have I3r -I- 1 unknown variables (fields) that are de¬ 
scribed with the quantities |na, [/^, Jg, T, Wg, q((, | 

(see Ref. 0). Such an expansion reads 

fa = [1 + Afpata + A'f'PafaPau] , (36) 

where fjf^ is the local equilibrium distribution function 
(Jiittner distribution) described by Eq. ((B]). In Eq. (I5B1) 
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the unknown tensorial coefficients are calcu¬ 

lated by solving an algebraic system constructed with the 
help of the definitions of the particle four-flow and the 
energy-momentum tensor . The details of such cal- 

I 


culation are long and unnecessary to do them here, they 
can be consulted in Refs. [1, [l^. As a result of such 
development the distribution function fa will depend on 
linear terms of the thermodynamic fluxes, namely. 


f _ f( 0 ) J 1 _ I Pa 


CaGa- 


u^p: 

kT 


Pa{fii') Ca 
2pa rna\^ 


Li V 

-PaPa 


+ 


Wg d In Ca 
Pa '91ncS 


a ' 'a 

3(c“ + K/T) U^pif clC + 3(c“ - h^/fcT^) 


fc2T2 


kT 


(37) 


Here we have introduced the abbreviation Gg = 
ATa (Ca) IK 2 (Ca) and the partial specific heats per par¬ 
ticle z% = k [C + 5GaCa - GlC - 1 ) and c“ = c“ -P fc at 
constant volume and pressure, respectively. 

Then, following the combined Chapman-Enskog-Grad 
method [^, the expansion (1571) must be compatible with 
the truncated Chapman-Enskog series up to first order, 
that is fa = {1 + 4 ) 0 ) where 4>a is the first correction 


to the distribution function . 

Now we can proceed to linearize the Boltzmann equa¬ 
tion as follows, we substitute fa = f^'’ (1 -I- 4a) in the 
left hand side of the Boltzmann equation ([5]) and keep 
the linear terms. This process is technically the same 
of that developed in the Chapman-Enskog method. We 
also use the so-called functional hypothesis, namely, fa = 
fa{xf^,pf\na,Uf^,T), leading to 


P)) 


df 


(0) 




-T! 




df, 


(0) 


dp>^ 




P'i ^ 

Ha dx>^ T 


1 ~ CaGa -b 


PaU\ 

kT 


C PaPaPjS^jdkl dgi X 

2kT U'^Par 


dr r ^ 


_ dxf^ 

PaPaPi^)j 1 
U'^Par J ' 


(38) 


On the other hand, we substitute the Grad function Boltzmann equation ([5|), and keep only the linear terms. 
Eq. (I37|) in the collisional term (right hand side) of the This process yields 


E [ (fan - fah)F,aaab = ” E Kl ^ + ^a, K] ^ 


-T 


ah 


kT 


Pt (. r ^^Pb 

T 0 ‘ ^ 


—Tab 


Pb/i q- 

TPb 


4 c.G.-bA 



Pa/i 

T , 

Ch 

jj 

TPa 

-f-'ab 



Pb{fii' 

2P6 



Pa+u) ^ 

91nCt- 

( UKJi.PbP'/, 

3(c^ + hb/r)c/xY 

ZUb 

1 PaPa 

[maK J 

„ Xah 

2Pa 

91nc(( 

y k^T^ 

K kT J 

Pb 


-I, 


ab 


5 In Ca / 3(c“ + K/T) U^p>) 


9 In c“ 


fc2T2 


kT 


Pa 


(39) 


Here we have introduced the collision operators 


^ab 


[Pa] = J n°^fb'’ {p'a - Pa) FabCJabd^y/^ 


d^Pb 

PbO 


(40) 

for any function that depends on the momentum four- 


vector ipa{Pa)- Note that Eq. (HOI) imply that we can 
write for an arbitrary function fjbiPj/) 


IpbTab [Pa] ^ = [ PaTab Kb] ^ 

PaO J PaO 


(41) 
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thanks to the symmetry properties of the collision oper¬ 
ator. 


By collecting the above information the linearized 
Boltzmann equation in the combined Chapman-Enskog- 
Grad method becomes 


f(0) I Pa 9na 


K 

T 


1 ~ CaGa + 


kT 


dT 1 i dUj PaP'aPi^rjhl dffl P^PaPiSi: 

dx>^ 2kT U^par dr r kT^^ U^par 


= - ^ K] ^ K] ^ - lab 
6=1 I Pb Pa 


Pt (p r ^•^Pb 


^bfj, q- 


■ah 


^ CgGg- 


kT 


—Tab 


-I, 


■ah 


a In Ca / t/^C/,pK 3(c“ + K/T) U^pii 


i91n c“ 


A:2t2 


kT 


^g/i 

Wa 



Pb{pu) ^ 


Pa(//i/) 

T , 

91nC6 

( U^U.p>^pl 

3(c^ + h6/r)[/xY 

ZUh 

mbhb^‘’^‘’ 

2P6 

1 PaPa 

[maha J 

2Pa 

-^ah 

u 

1 

y k'^T'^ 

4 kT J 

Pb 


^ , (42) 


due to (l38|l and (|3^ . 

In the next sections we will use (H^ in order to deter¬ 
mine the constitutive equations for the diffusion fluxes 
J((, heat flux q^, non-equilibrium pressure w and pres¬ 
sure deviator tensor pG'^). 


IV. FICK AND FOURIER LAWS 

Now we will obtain a system of linear equations for 
the determination of the the diffusion fluxes J)) and the 
heat flux of the mixture q'^. The solution of such a sys¬ 
tem will represent the form of the linear fluxes in terms 
of the thermodynamic forces. The integral functions for 
the transport coefficients and therefore the Onsager reci¬ 
procity relations will be analyzed in the next subsection. 

To obtain the first one of the looked set of equations, 
we multiply Eq. by cA((p((/nQ and integrate over 




d^'pa 

PaO 


The integrals used for this process can be 


consulted in the Appendix. The resulting equation is 




> am* 


U‘ 


1 


9$ 


.m _ 

dx’' 1 — 4)2/4(34 ^2;* 

r 

= ^ (Aab^, - AabP^) (43) 
6=1 

where = A^’^di, is the gradient operator and $ = 
— is the gravitational potential. In Eq. (l43l) we have 
introduced the matrices Aab and Tab- We can split Aab 
for different indices {a,b} 


Aab — 




PapTab [l)6z/] "x/ "P' 


SnantkT 
and for equal indices {a, b = a}, 
cA'" 


d^Pa 

PaO 


a ^ b, 
(44) 


T PapTaa \p au\ 




d^Pa 

PaO 


(45) 


The matrix Tab introduced in Eq. (1431) is written by 
doing the same splitting, for unlike indices {a, 6} we have 


Tab = — 


cAf^- 


SnaObkT^ 


PapTab 


Cb UrPl , 

-T G6- a Pbv 

c“ V rubT 




d?Pa 

PaO 


a^b, (46) 


and for like indices {a, b = a}, 

cA^^’^ I 

^ 1 

6=1 ■ 


3n2fcT2 


Pa^^a 



^ (c - I p 

.r, \ ___ _o I Pai' 


rUaC^ 


I ^^.(47) 


PaO 


Next we look for a second equation which is indepen¬ 
dent from Eq. (1431) . Hence, we multiply the linearized 

Boltzmann equation ^ by A^-^0^ (Ca - K 

and integrate over x/—5^^^, for this long process we use 
also the integrals that appear in the Appendix. The re¬ 
sult becomes 


If T . 

^2 


T 


W' 


1 


94) 


m _ 

dx’' 1 — 4)2/4c^ 9a;* 

r 

= Y,{TbaK -'^abC^^b) m 

6=1 

where another matrix T-Lab is defined. As with the others 
operators, we split T-Lab in the part for unlike indices {a, b} 


Tiab = — 


cA'*** 

3nan6fcT3 




TOaC* 


A — - 

'A^aa — 


3n2fcr 


^ ^ / Pafi^ah [Pai^] 

^^ah 

\ Cb u,pi \ 1 

,, G6 2 IPb'' 

L 6=1 


[c® \ rubT J \ 




d^Pa 

PaO 


a ^ 6,(49) 
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and the corresponding for same indices 


^'^ab 


'Haa = — 


^ G.- 


xl„ 



Hence, we have obtained the desired system of alge¬ 
braic equations, namely (l43l) and (l48l) which are an inde¬ 
pendent set of linear equations for the determination of 
the diffusion and heat fluxes. 


force as 

- 1 —1 r 

= (52) 

6=1 6=1 

Above we have considered the constraint J2a=i ^ 
which implies that there exist only r — 1 linearly inde¬ 
pendent diffusion fluxes. 

Next we recall that the chemical potential of species 
a is defined through the Gibbs function per particle 
(/ia = ^a — TSa + Pa/^a)- So that, the following important 
relationship holds for its gradients 

( y) = 

Therefore, the substitution of Eq. (l53l) into (l43l) yields 


A. Onsager reciprocity relations 


- { y ) ~ ^ ^ 

6=1 


In this section we will show that the Onsager reci¬ 
procity relations hold for the system under consideration. 
The idea is to verify if the matrix associated with the 
diffusion coefficients is symmetric and therefore the so- 
called cross effects are equal as it is described from one of 
the hyp othesis of the Linear Irreversible Thermodynam¬ 
ics [1^. One cross effect for our system is the contribu¬ 
tion to diffusion due to the temperature gradient, this is 
often called “Soret” effect. The other cross effect is the 
contribution to the heat flux due to the chemical poten¬ 
tial gradient or a concentration gradient, when it is due 
to the last, it is called “Dufour” effect. This demonstra¬ 
tion is general in the sense that no interaction between 
the particles is established, but of course, the microscopic 
reversibility principle is called for the collisional term of 
the Boltzmann equation ([5]). 

Let us now write the thermodynamic forces in order to 
identify clearly the Soret and Dufour effects in terms of 
the temperature and chemical potential gradients. 

First we define a generalized thermal force as 


dUi I 5$ ■ 

dx'' 1 — $^/4c^ dx'^ ’ 

(51) 

where the first term contains a temperature gradient 
while the second one - whose nature is strictly rela¬ 
tivistic due to the factor Tjb? ~ is proportional to the 
four-acceleration and the gravitational potential gradi¬ 
ent. The term due to the four-acceleration was proposed 
by Eckart [l^ while the one due to the gravitational 
potential gradient by Tolman [mil- If we think in 
a relativistic gas in equilibrium, we can conjecture the 
following two aspects: (i) in the absence of a gravita¬ 
tional potential gradient, the temperature gradient must 
be counterbalanced by an acceleration and (ii) in the ab¬ 
sence of an acceleration, the temperature gradient must 
be counterbalanced by a gravitational potential gradient. 
Now, Eq. (H51) can be written in terms of the thermal 


V^T = V'^T - —A'"* 
(? 


Moreover, by considering that there exist r — 1 indepen¬ 
dent diffusion fluxes, we we can take the rth component 
of Eq. (IMl) and subtract it from (IMl) itself, yielding 

— ~ Mr ^ ~ ^r J- _ 

r—1 

= ^ ^ -^ar “1” J^rr) 

b^l 

r 

b^l 

Now we can proceed to solve the system of linear equa¬ 
tions formed by Eqs. (|52|l and (|5^ . First we solve Eq. 
(f52l) for q^, yielding 

d^l ^ ^ 

r r—1 

+T.T.i'^~")cA^bd-j^rd)r,, (56) 

d—1 b—1 

where is the inverse matrix of T-Lcd so that 

{n-^),d'^da = Sea is the identity matrix. Then we insert 
Eq. (1^ into Eq. (1551) and solve for J((, 


r— 1 

Tj2Kb^^ 

6=1 


fib f^r 

f 


T> 

—v^r. 

T 


(57) 


Here we identify the above equation as the generalized 
Fick law, where the coefficients and Va are related 
with the diffusion and thermal-diffusion (Soret) effects, 
respectively. The inverse of the diffusion matrix reads 

I'D = A-ab ~ Arb ~ Aar + Arr 

V / ab 

r r 

- E E i'^~')cd (^bd - J^rd ), (58) 

c—1 d—1 


















while the thermal-diffusion coefficients are given by 

r—1 I" r r 

\^b-K+YY i^bc - Trc) (H’') 


cd 


b^l 


c—1 d—1 


(59) 

Now we have to obtain the total heat flux as a function 
of the temperature and chemical potential gradients. For 
this end, we rewrite the total heat four-flux (1^51) as 


a=l 


r— 1 

E 

a=l 




(60) 


and substitute in it the expressions found for q)) and J)), 
i.e. Eqs. dSH) and (EH). Hence it follows the Fourier law 


A' 


r —1 


= --v^r-r^2?;v^ 


a=l 


/^a Mr 

f 


(61) 


where we have introduced the thermal conductivity co¬ 
efficient 


r r r—1 

a—1b—1 b—1 


hh - K 


EE(^ ^)ac^^bc J^rc) 


a—1c—1 

and the diffusion-thermal coefficient 


( 62 ) 


r—1 




6=1 


hh - hr 


c—1 d—1 

(63) 

Ultimately, we make a close inspection of the matrices 
Aab, J^ab and T-Lab which are given as functions of the col¬ 
lision operators lab- From (l44l) . (l46l) and (|4^ . we may in¬ 
fer that only ylab and T-Lab are symmetric matrices, while 
J-ab is non-symmetric. Hence we may conclude from (1551) 
that the matrix related with the diffusion coefficients are 
symmetric, i.e., Moreover, for the coefficients 

of cross effects - namely the Soret Va and Dufour - 
we note from the symmetry of Tlab and 'Dab that (1591) and 
(l63l) are equivalent, so that Da = D'^^. The relationships 
D'^f^ = D'^^ and Da = D'^ imply a demonstration of the 
validity of the Onsager reciprocity relations. Note that as 
in another’s demonstrations [2^, it appears to be an ulti¬ 
mate macroscopic effect that can only be proved because 
of the symmetries that belong to the collisional term of 
the Boltzmann equation dS]) given from the H-theorem 
i.e. microscopic reversibility principle. 


B. Thermal and diffusion forces 

It is usual in the theory of fluid mixtures to express the 
diffusion fluxes and the heat flux of the mixture in terms 
of the generalized thermal and diffusion forces. The ther¬ 
mal force for a relativistic fluid was introduced in the last 


section (see (ISTI) !. On the other hand, we follow Q and 
define the generalized diffusion force of species a as 


d^ = V'^Xa + (Xa-l)V'^lnp 

nghg- nh - f _ 1 

pc^ dx'^ 1 — <1)2 dxT 


( 64 ) 


where Xg = p^/p = Hg/n is the concentration of species a. 
We can identify four contributions to the generalized dif¬ 
fusion force: a concentration gradient, a pressure gradi¬ 
ent, a term proportional to the four-acceleration and the 
gradient of the gravitational potential. Here it is impor¬ 
tant to emphasize that contrary to what happens with 
the thermal force, the terms with the four-acceleration 
and gradient of gravitational potential are not of a 
strictly relativistic nature. Indeed, (nghg — nh)/c^p = 
(ngmgGg - nbm{,Gb)/p and Gb 1 for (b » 1- 

This equation has a very important feature because it 
represents the generalization of the diffusion force origi¬ 
nally written for the non-relativistic case [U [3 ■ More 
discussions about this point can be found in [7|. Here 
as in the non-relativistic case, exist only r—1 linearly 
independent generalized diffusion forces due to the rela¬ 
tionship Y^a=l 6a = 0. 

Now we can proceed to express the vectorial fluxes in 
terms of the generalized thermal and diffusive forces. To 
do sm we use the momentum density balance equation 
(see 0) 


dp nh 
dx^ (? 


U‘ 


. dUj 1 d<^> 

dx'^ 1 — < 1 ) 2 / 40 "^ dx'^ 


= 0 , 


(65) 


and the gradient of the chemical potential written as 


- TV'" 


f^a Mu 
T 


= - — (V'^Xg -h Xg V'" In p) 


+ -^ (V'^x, + x.V'" In p) + ^ V'^T. (66) 

Dr 1 

After some rearrangements, the expressions for (I57p and 
dM]) become 


r _(L 

= + (67) 

6=1 

q'" = + (68) 

a—1 

In these representations for the generalized thermal and 
diffusion forces, the transport coefficients read: 


r— 1 


c=l 
r — 1 r r 


rib 




p. = - U W (x“‘) 


cd ' 


b—1 c—1 d—1 
r r 


a = -EE 


a=l6=1 


i^Yab 


(69) 

(70) 
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r— 1 

c=l 


V' 


r — 1 




(71) 

(72) 


At this point it is worth pausing to make two com¬ 
ments. Firstly, by looking the expression for the general¬ 
ized diffusion force Eq. (IMl) we note that it depends on: 
(i) A concentration gradient, that tends to reduce the 
non-homogeneity of the mixture; (ii) a pressure gradient, 
where heavy particles tend to diffuse to places with high 
pressures, e.g. in centrifuges; (iii) an acceleration, which 
acts on different masses and (iv) a gravitational potential 
gradient. Secondly, let us suppose a mixture in which the 
generalized thermal force vanishes, the pressure is con¬ 
stant and there is no acceleration. We can think also that 
there is no diffusive flux, implying a pseudo-equilibrium 
state. It is very interesting such a situation because of 
its physical implications, that is, the gradient of concen¬ 
tration has to be counterbalanced by the gravitational 
potential gradient. 

To complete this section, we point out that the thermal 
conductivity coefficient A in a mixture is defined as the 
ratio of the heat flux to the temperature gradient. This 
when there is no diffusion i.e., when = 0. From (1561) 
and dini), we have 


^ = EEMa.- (73) 

a—1 6^1 

Furthermore, in the absence of a temperature gradient 
the constitutive equation for the diffusion fluxes (1571) are 
proportional only to the generalized diffusion forces and 
(Dab is identified as the matrix of the diffusion coefficients. 


V. ENTROPY FLUX OF THE MIXTURE 


In this section we will show that the entropy four-flow 
for the system under consideration takes the form as pre¬ 
dicted by the Linear Irreversible Thermodynamics. Ac¬ 
cording to (1271) and dsniia the entropy flux of the mixture 
is given by 


= -kcAi^Y, 

a—1 


Pa/aln(ba/a)\/^ 


d^Pa 

PaO 


(74) 


We substitute the Grad distribution function dSZl) 
into the above expression and linearize in the fluxes 
J^, . After integration the entropy flux takes 

the form 


a=l a—1 

_ 



(75) 

a=l 


The first equality above shows that the entropy flux of 
the mixture is a sum of two terms: one refers to the 
sum of all partial heat fluxes divided by the temperature 
and the other is a sum of the transport due to diffusion 
of the partial entropies per particle. The second equal¬ 
ity is well-known from non-relativistic Linear Irreversible 
Thermodynamics , and is connected with the trans¬ 
port of the chemical potentials driven by diffusion. 

We can also express the entropy flux of the mixture in 
terms of the thermal and diffusion generalized forces by 
substituting the representations (1571) and (1551) into (1751) . 
yielding 

= (76) 

a—1 

Here the scalar coefficients C and Ca read 


r — 1 

C = X-'^{pa- Pr) Va, 

a—1 

(77) 

r — 1 

Ca=T>'^-'^ {fib - fir) i^ba- 
6=1 

(78) 


It is clear from the definitions of the thermal (15T|) and dif¬ 
fusive ((551) forces that the entropy flux of the mixture (1751) 
depends on the temperature, concentration and pressure 
gradients as well as on the acceleration and gravitational 
potential gradient. 


VI. NAVIER-STOKES LAW 

In this section we will calculate the constitutive equa¬ 
tions for a relativistic Newtonian fluid, in other words 
the Navier-Stokes law. This law is usually separated in 
two equations. The first one is for the non-equilibrium 
pressures and it is associated with the bulk viscosity. The 
second one is for the pressure deviator tensor and it is 
associated with the shear viscosity. 

Let us start with the constitutive equation for the par¬ 
tial non-equilibrium pressures vua of species a, it is ob¬ 
tained as follows. We multiply (l42)) by A^rPaPa 

_ ,3 

integrate over \/—g-^■ For this purpose, we use the in¬ 
tegrals from the Appendix. We also eliminate the deriva¬ 
tive projections with the help of the partial particle 

number density and energy per particle balance equa¬ 
tions. Such balance equations correspond to an Eulerian 
fluid, where non-equilibrium quantities J((, 
vanish, that is 

= 0, (79) 

n,c“C/^9^r+p,V'^17^ = 0. (80) 

The result of this process becomes a system of equations 
for Wb and it reads 

r 

V =Y,HabWb. (81) 
6=1 


p„fcr i91nCa 
d In eg 
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Here we have introduced the matrix 'R,ab^ which is defined 
for different indices {a,b} as: 


l^n.h — 



And in a similar fashion of the previous sections, we write 
this matrix for equal indices {a, 6 = a} as 


T^aa = 


C^Pa 


/ PaPa^a 




ainCa / UrPl 


d In c“ 


kT 


3(c“ + hjr)\p- 


+ / p>:ia 


eg J kT 

3(c° + K/T) \ p- 
eg / kT 


d In Ca ( UtPI 
d In eg 


kT 


>V^^. (83) 

PaO 


The solution of the linear system of equations (ISTl) for 
the partial non-equilibrium pressures vJa is given by 


'^a 


r 




PfafcT c)lnCb 
d In eg 




(84) 


where (TZ ^)^j^ denotes the inverse of the matrix TZab- 
The constitutive equation for the non-equilibrium pres¬ 
sure of the mixture is obtained from the sum of ((5^ over 
all constituents according to ([231)3. Hence it follows 

et7 = (85) 


where the bulk viscosity coefficient of the mixture reads 


^ = E ")ab 

a,b—l 


PfcfcT dln(b 
d In eg 


( 86 ) 


For the second equation that conforms the Navier- 
Stokes law, which is the pressure deviator constitutive 
one, we proceed in an analogous manner. We take the 
product of with — A^r A'^‘'/3]p^p^ and inte- 

,3 

grate over This process leads to the following 

linear system of equations for the partial pressure devia¬ 
tor tensors : 


2v(m;7-> (87) 

6=1 


In this last equation we have introduced the following ab¬ 
breviation for the symmetric and traceless four-velocity 
gradient 




A^Ag -p AgAg A'^'^Ao 


d'^U^.{ 88 ) 


Equation (l87l) also includes the definition of the matrix 
ICab which reads 


ICab = — - 


lOp^haPb 


ICnn. = -- 


C^A„/„A 




lOPahaPa 


PZPl^ab 

, - d^Pa 

PaO 

E / P^’ofJ-ab 
6=1 ■! 


Cb u u 
a (89) 

^a n 1/ 

- irP^Pa 


1 PaPldlaa 

n y 

. PaPa 

1 



i 


(90) 


PaO 

From the solution of the linear system of equations (l87l) 

for and from the relationship = X]6=i P6^^^ 
follows the constitutive equation for the pressure deviator 
tensor of the mixture: 

(91) 

Here the shear viscosity coefficient of the mixture is given 
by 




(92) 


a,b—l 


Equations (l85l) and (IMI) are the constitutive equa¬ 
tions of a relativistic Newtonian fluid, also known as the 
Navier-Stokes constitutive equations. 


VII. CONCLUSIONS 

In this work we have studied a mixture of r species of 
relativistic gases in the presence of gravitational fields. 
The curvature of the space-time was introduced by in¬ 
corporating the Christoffel symbols to the Boltzmann 
equation. We used the Schwarzschild metric written in 
isotropic coordinates. A linearized Boltzmann equation 
was obtained by following a methodology which combines 
the features of the Chapman-Enskog and Grad methods. 

By applying the Chapman-Enskog-Grad combined 
method to the Boltzmann equation we obtained a linear 
expression (Eq. (|42l) l which was used for the determi¬ 
nation of the thermodynamic fluxes as functions of the 
thermodynamic forces. The Navier-Stokes law was de¬ 
rived as well as the generalized of Fourier and Fick laws. 

The proof of the validity of the Onsager reciprocity 
relations was possible lastly because of the symmetries 
of the collisional term of the Boltzmann equation. These 
symmetries are those associated with the H-theorem and 
the microscopic reversibility principle. This reinforces 
the idea that the Onsager reciprocity relations are the 
macroscopic manifestation of the microscopic symmetries 
of the trajectories of the particles that conform the gas. 

We have introduced the thermal force 

V'"T = V^T- — A^* 

(93) 


dUi 1 

dx'^ 1 — /4:C^ dx^ 
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which is a very eloquent result. Indeed, Eq. (I93p turns 
to be just the gradient of the temperature in the non- 
relativistic limit i.e. V^T because the factor T /c^ of the 
second term is of relativistic order. The inclusion of the 
acceleration term into the thermal force was proposed 
by Eckart [H while the one relating the gravitational 
potential gradient by Tolman MM- Here these terms 
appear as a natural consequence of the solution of the 
relativistic Boltzmann equation in gravitational fields. 

On the other hand we have identified the generalized 
diffusion force with 


d^ = V'^Xa + (xa-l)V'^lnp 


nghg - 
pc^ 


"'If 




1 — $2/40"^ dx^ 


■ (94) 


This is a new and interesting result because the third 
term ~ which is related with the four-acceleration and 
the gradient of the gravitational potential - does not go 
to zero in the non-relativistic limiting case as the thermal 
force. As it was pointed out in the work Q , the diffusion 
force that came out from a non-relativistic kinetic the¬ 
ory [l3, [3 has a similar expression to (IMl) . It depends 
on the concentration and pressure gradients, but it has 
a term depending on the forces which act on the parti¬ 
cle of different species, such a term vanishes when only 
gravitational forces are acting on the particles. 

Another result obtained is the entropy flux of the rel¬ 
ativistic mixture through the use of Grad’s distribution 
function, which has a similar expression as the one of non- 
relativistic Linear Irreversible Thermodynamics M- Its 
constitutive equation was written in terms of the gen¬ 
eralized thermal and diffusion forces, so that it depends 
also on the acceleration and on the gravitational poten¬ 
tial gradient. 

Here is the place to discuss two additional issues. The 
first one is the validity of the Onsager reciprocity rela¬ 
tions for the case of a relativistic quantum gas. In such 
a case, the local equilibrium distribution (which in this 
work is given by Eq. ©) would take a form of the Eermi- 
Dirac and Bose-Einstein distributions for a fermionic and 


bosonic gas, respectively. Quantum relativistic gases can 
be described by the relativistic Uehling-Uhlenbeck equa¬ 
tion (see e. g. Q)- As we have pointed out, the validity of 
Onsager’s reciprocity relations are deeply associated with 
the symmetries that belong to the H-theorem. In the 
present work, those symmetries are implied in Eq. m- 
Then, to show the validity of the Onsager reciprocity re¬ 
lations for a quantum system we need the validity of the 
H-theorem, fact that has been shown in the literature in 
Refs, and raising the possibility to explore that 
issue. The second topic is related with Tolman’s law, 
which has been derived in [I^ and is valid for all static 
spherical symmetrical line element. In the present work 
we have used the Schwarzschild metric, which accord¬ 
ing to Birkoff’s theorem is the most general spherically 
symmetrical non-rotating and uncharged source of the 
gravitational field. 

As a final comment we call attention to the fact 
that for the determination of all the transport coef¬ 
ficients, we have to specify the interaction potential 
of the relativistic particles and evaluate the matri¬ 
ces {TZab,}Cab,Aab,J^ab,'Hab}- This represents work in 
progress and will be published in the future. 
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APPENDICES 
Table of integrals 

For the purposes of this work, it is convenient to do the 
following unique decomposition of the integral operators: 


Kb = J p'^iab m ^ {Kb - lab) i (4\ - K) 5^^", 


PaO 

.d 

PaO 


Kr = I pWab K] ^ {Kb - K) 


+Y^ilb-ilb) 


Kr = j p^aiab [KK] rg^ = ( 244 \ -124^ + K) + ^(/J, - 2K 


-Kb)K'^U'^U- + g^^u^un + ^ ( 16 /®, - 3 / 1 , - 12 /®,) 

^ ( 3 j 7 ^ _ + 2 /®,) + g^^^g^ + ^ {K - K + Kb) 


(95) 


(96) 


( 97 ) 
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where ... /Jf, are given by 


_ Ufj.Ui' 

^ab 2 ^ab ’ 


t 2 _ tS _ UfJ.^vU<T yHV(T 

^ab dtiv^ab ’ ^ab ^3 


jb _ UyJJyUfjUr 

^ab - -4 


^ab ’ 


j-6 QurU^U(j 


4 _ Q^KjUy 


r = 

^ab 




ab ’ 


-y^uar t 7 _ ^ ^ yfivar 

^ab ’ ^ab “ 9^<Tyi^T^ab 


(98) 

(99) 


I- 

Here we list a table of integrals that are used in the 
previous sections. 


f e kT^ ^ = 47rTOfcriiri(C), [p^e —E.=4TYm^kTK2{C)U^, 

J Po J Po 

/ I A r 

-= -47r(TOfcr)^ K2{C,)g^^ - C^siC) — 5 - 

Po L 


p^p'^p'^e kT^>-P = —47rm^(fcr)^ 
Po 




p'^p p p e 


Po 


3 

Km 




p^p'^p'^p'^p’^e = 47rm^(fcT)^ 


Po 


3 ^ ^ *v,. g^2 


+C"^6(C) 


UPU'^U'^U'^U^ 


e Tk'U’'px (pp 

U^Pr Po 


/ g l^'J' U Px fQ 

pP —^-7^-= 47rm^ 


^i(C) 




e ik'U px (]3p Airm^kT 

U'^Pr Po 3 


U^Pr Po 

[K2(C)-C(^i(C)-Kii(C))]p'^^ 




-l[4/f,«) - c(A',(c) - Kii(c))]t;'-c;"|, 

- [cKm) +2a,(o]11^|. 

J U^Pr Po ( 3 ■' J 

d^p [ 3^3(0 - Cj^2(C) + CiKijO - Kii(C)] 


pf^p'^p'^p^ 
1 


U^Pe Po 


15 


1 


9"" 9 


[18^3(0 - CK^iO + e {Ki{0 - Kii(C))] 


+- 


6c^ 

3 


48i^3(C) + 4Ci^2(C) + C" (^i(C) - Kii(C)) 


U^V'U'^U'' 


pf^p'^p’^p^p^ 


e kT^ PA (]3p AirmPc^ ( ii^3(C) 

U^Pe Po I 15 




9 9 


«) _ 


10c2 


[ 8 X 3 ( 0 +CK2{0]9^^^U''U'^U^ 


■[C^KsiO + 12CK2(0 +80K3(0]Uf^U'^U'^U^U^\. 


( 100 ) 

( 101 ) 

(102) 

(103) 

(104) 

(105) 

(106) 
(107) 


(108) 


(109) 


Above the parenthesis around N indexes indicate a sum over all permutations of these indexes divided by N\. 
Furthermore, Ki„(C) denotes the integral 


Ki„(C) = / 

Jo 


00 g —C cosh t 

0 cosh "t 

I- 


dt. 


( 110 ) 


Cristoffel symbols for the Schwarzschild isotropic 
metric 
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pU _ r, pU _ n 

^ 00 ~ '-'j ^ ii ~ 


pO _ 1 dgo{r) ^ 

2go(»') dr r 


rij=0 {i^j^k), 
1 


ro, = 0, 


pz - dgo{r) 

°° 2 gi{r) dr r 


J dgijr) ^ 

2 gi{r) dr r 


r-. = 

13 

TL = - 


1 dgi{r) x^ 
2 gi{r) dr r 




( 111 ) 

( 112 ) 


The underlined indices above are not summed and 

dgo(r) 2GM {l - §£) 

dr c^r"^ ( 2 ^ _|_ 


dgijr) 2GM / GM^ 

dr c^r^ \ 2c^r ) 


(113) 


[1] F. Jiittner, Das Maxwellsche Gesetz der 
Geschwindigkeitsverteilung in der Relativtheorie, 
Ann. Physik und Ghemie 34 (1911) 856. 

[2] S. R. de Groot, W. A. van Leeuwen, Gh. G. van Weert, 
Relativistic Kinetic Theory, North-Holland, Amsterdam, 
1980. 

[3] G. Gercignani, G.M. Kremer, The Relativistic Boltzmann 
Equation: Theory and Applications, Birkhauser, Basel, 
2002 . 

[4] A. Sandoval-Villalbazo, A. L. Garcia-Perciante, D. Brun- 
Battistini, Tolman’s Law in Linear Irreversible Thermo¬ 
dynamics: a Kinetic Theory Approach, Phys. Rev. D 86 
(2012) 084015 

[5] M. Smerlak, Diffusion in Gurved Spacetimes, New Jour¬ 
nal of Physics 14 (2012) 023019. 

[6] G. M. Kremer, Relativistic Gas in a Schwarzschild Met¬ 
ric, J. Stat. Mech. (2013) P04016, E05001, ElOOOl. 

[7] G. M. Kremer, Diffusion of Relativistic Gas Mixtures in 
Gravitational Fields, Physica A 393 (2014) 76. 

[8] W. Zimdahl and G. M. Kremer, Temperature Oscil¬ 
lations of a Gas in Gircular Geodesic Motion in the 
Schwarzschild Field, Phys. Rev. D 91 (2015) 024003. 

[9] A. G. Bezerra Jr., S. Reinecke, G. M. Kremer, 
A Combined Chapman-Enskog and Grad Method. 1. 
Monoatomic Gases and Mixtures, Continuum Mech. 
Thermodyn. 6 (1994) 149. 

[10] G. M. Kremer, On the Kinetic Theory of Relativistic 
Gases, Continuum Mech. Thermodyn. 9, (1997), 13. 

[11] H. Grad, On the Kinetic Theory of Rarefied Gases, Com- 
mun. Pure Appl. Math. 2 (1949) 331. 

[12] C. Eckart, The Thermodynamics of Irreversible Pro¬ 
cesses, III. Relativistic Theory of a Simple Fluid, Phys. 


Rev. 58 (1940) 919. 

[13] R. C. Tolman, On the Weight of Heat and Thermal Equi¬ 
librium in General Relativity, Phys. Rev. 35 (1930) 904. 

[14] R. C. Tolman and P. Ehrenfest, Temperature Equilib¬ 
rium in a Static Gravitational Field, Phys. Rev. 36 (1930) 
1791. 

[15] C. W. Misner, K. S. Thorne, J. A. Wheeler, Gravitation, 
W. H. Freeman and Company, San Francisco, 1973. 

[16] S. R. de Groot, P. Mazur, Non-equilibrium Thermody¬ 
namics, Dover Publications, New York, 1984. 

[17] S. Chapman and T. G. Cowling, The Mathematical The¬ 
ory of Non-Uniform Gases, 3rd. edition, Cambridge Uni¬ 
versity Press, Cambridge, 1970. 

[18] J. O. Hirschfelder, C. F. Curtiss, R. B. Bird, Molecular 
Theory of Gases and Liquids, John Wiley, New York, 
1964. 

[19] G. M. Kremer, W. Marques Jr., Grad’s Moment Method 
for Relativistic Gas Mixtures of Maxwellian Particles, 
Phys.Fluids 25 (2013) 017102. 

[20] V. Moratto, A. L. Garcfa-Perciante, L. S. Garcfa-Colm, 
Validity of the Onsager Relations in Relativistic Binary 
Mixtures, Phys. Rev. E 84 (2011) 021132. 

[21] E. C. G. Stueckelberg, G. Wanders, Thermodynamique 
en Relativite Generale, Helv. Phys. Acta 26 (1953) 417. 

[22] H. Hebenstreit, Balance Equations for a Relativistic 
Plasma, 1. Differential Term, Physica A 117 (1983) 631. 

[23] G. M. Kremer, C. H. Patsko, Relativistic Ionized Gases: 
Ohm and Fourier Laws from Anderson and Witting 
Model Equation, Physica A 322 (2003) 329 . 

[24] C. Gercignani, G. M. Kremer, On relativistic collisional 
invariants, J. Stat. Phys. 96 (1999) 439. 

[25] C. Gercignani, G. M. Kremer, Trend to equilibrium of a 
degenerate relativistic gas, J. Stat. Phys. 98 (2000) 441. 













